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Abstract 

This is an exposition for mathematicians of some unsolved prob- 
lems arising in control theory of linear time- independent systems. 

The earliest automatic control devices that I know are described in the 
book of Hero of Alexandria "Pneumatica" , see Fig. 1. In the modern times 
these devices are omnipresent (almost every home appliance contains at least 
one, a car has several, an airplane or a guided missile has many; an ingenious 
mechanical steering device of a sailboat permits you to sleep and to dine 
during your voyage, while it keeps prescribed direction with respect to the 
wind; one can add many other examples). 

The mathematical theory of these devices begins, as far as I know, with 
George Biddell Airy (of the Airy function), Astronomer Royal, who inves- 
tigated mathematically stabilization of the clockwork mechanism directing 
his equatorial^ The stability condition that "all poles must be in the left 
halfplane" was explicitly stated for the first time by J. C. Maxwell [20]. Par- 
allel research was done in Eastern Europe by Aurel Stodola (1894) and Ivan 
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1 Equatorial is a device that continuously adjusts a telescope direction to compensate for 
the diurnal rotation of the Earth. One of the most complicated modern control systems 
directs the Hubble telescope. It has to keep the direction of the telescope with high 
accuracy and for long time. 
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Vyshnegradsky (1877) who pioneered the use of complex function theory, 
anticipating the work of Nyquist (1932), see for example, [TB] . 

Most of the XIX century research in the area was related to stabilizing 
the system which consists of a steam engine controlled by the governor. 

These investigations led to the famous criteria in terms of coefficients of a 
polynomial for all its roots to belong to the left half-plane, (E. Routh, 1877, 
A. Hurwitz, 1895), see [TU]. 




Fig. 1 A XIX century illustration made according to the description in the 

book of Hero Pneumatica. 

A linear system of the simplest kind is described by 3 real matrices: 
(A, B, C) of sizes n x n, n x m and p x n respectively. We have vectors 
depending on time: the inner state x(t) with values in R n , the input u(t) 
with values in R m and the output y(t) with values in R p . These are related 
in the following way: 

x' = Ax + Bu, 
y = Cx. 

We will only consider the case m = p = 1 (so called single input - single 
output systems). 
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Taking Laplace transforms, and assuming that x(0) = u(0) = y(0) = 0, 
we obtain zX(z) = AX(z) + BU(z), Y{z) = CX(z), so 

Y{z) = C{zl - A)- 1 BU(z) = p{z)U{z). (1) 

The rational function p(z) = C(zl — A)~ 1 B is called the transfer function. 
It is real and p(oo) = 0. Rational functions satisfying p(oo) = are called 
proper. For every proper rational p function there exists a triple (A, B,C) so 
that p(z) = C(zl - A)- 1 B. 

The correspondence between triples of matrices and rational functions is 
not trivial, not bijective, and there is a large literature on recovery of A, B, C 
from the transfer function (realization theory). But all essential properties of 
the system are encoded in the transfer function and here we identify a linear 
system with its transfer function. 

Improper transfer functions are equally important, they arise from more 
general systems of differential equations with constant coefficients; I don't 
go into detail, but the primary object in this paper will be an arbitrary real 
rational function; we call it a transfer function. It completely describes a 
linear system. 

A system is called stable if the transfer function has no poles in the open 
right half-plane H. The poles of the transfer function are nothing but the 
eigenvalues of the matrix A of the system. 

For a given (maybe unstable) system, one may wish to stabilize it by 
attaching a feedback controller. A controller is a linear system of the same 
kind; it is described by another real rational transfer function c{z). Attaching 
a controller as in the third diagram in Fig. 2 means that we take the output 
of our original system, transform it by the controller, and then add to the 
input: 

Y = p(U + cY) = P U + pcY. 

We obtain a new system, which is called the closed loop system. By solving 
with respect to Y we get the closed loop transfer function: 



Cancellation between poles and zeros of c and p is possible here, but engineers 
naturally do not want to rely on such cancellation. So they give the following 
definition: 
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A controller c internally stabilizes p if 1 — cp has no zeros in the right 
half-plane H, the poles of c are disjoint from the zeros of p in H, and the 
zeros of c are disjoint from the poles of p in H. 

From now on by "stabilization" we mean "internal stabilization". One 
can easily show that internal stabilization is equivalent to the condition that 
all four transfer functions 

pc/(l-pc), c/(l-pc), p/(l-pc), 1/(1 -pc) 

are without poles in H. 

All these four transfer functions can be realized by attaching the feedback 
in various ways, as shown in Fig. 2 below. 





c 




p 












A 










Fig. 2 Configurations corresponding to 
pc/(l-pc), c/(l -pc), p/(l-pc), 1/(1 -pc). 
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Here is another elegant way to rewrite the internal stabilization condition: 
c internally stabilizes p iff c avoids 1/p in the sense that 

c(z) + l/p(z), z e H. (3) 

Thus the stabilization problem is: for a given rational function p find a 
rational function c so that (Tj|) holds. 

We obtain an equivalent problem when the right half-plane H is replaced 
by the unit disc D. From the point of view of system theory, the unit 
disc setting corresponds to discrete-time systems. Instead of a differential 
equation, we have a recurrence relation, 

x{n + 1) = Ax{n) + Bu{n), y{n) = Cx(n), 

and in place of the Laplace transform we use the generating function X(z) = 
X^oo x(n)z n . Then the transfer function p(z) becomes C{z~ 1 I — A)^ 1 B, 
exactly as in the case of continuous time, and the system is stable if p has 
no poles in the unit disc, which now is equivalent to saying that A has no 
eigenvalues whose absolute value is greater than 1. 

Stabilization of one system is always possible if one does not restrict the 
degree of c. 

Now we consider simultaneous stabilization of several systems by one 
controller. The problem has evident practical meaning: the system that we 
want to stabilize may work in several different regimes (think of the cool- 
ing/heating system in your home, which is usually controlled by a single 
devise, or an airplane during take of/landing/horizontal flight), mathemati- 
cally this means that we want a single controller to stabilize several systems. 

Consider the problem of stabilizing two systems p\ and p2- This is equiv- 
alent to stabilization of one system by a stable controller [3J. 

Stability of the controller is a desirable property by itself: if the system 
p suddenly stops working, we don't want the controller to destroy itself. 

A rational function avoiding two different rational functions pi in H al- 
ways exists, but it may have complex coefficients, even if p, are real. One 
usually needs a controller with real coefficients. There is an obvious topo- 
logical obstruction to the existence of real rational function without poles 
on R>o which avoids a given real rational function. In fact, this is the only 
obstruction: 

For a given real rational function p, there exists a real rational function 
c satisfying (T3|) and without poles in H if and only if p has even number of 
poles between every two adjacent zeros on R>o- 
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This neat statement is due to Youla, Bongiorno and Lu [22], and control 
theorists are very fond of it [3]. 

Now we consider simultaneous stabilization of three systems. 

Finding a function which avoids three given functions is an interesting 
problem which attracted attention of pure mathematicians who were unaware 
of its application to control theory. It seems that the problem was stated for 
the first time in [23], and a connection with an "interpolation problem" of 
the kind stated below in Theorem 1 was established. 

In [IT] this avoidance problem is considered for meromorphic functions 
which avoid given rational functions in an arbitrary given region. The author 
credits Volberg and Eremenko who stated the problem and obtained some 
partial results. Apparently they were motivated by the "Lambda-lemma" 
and holomorphic motions which were discovered about that time [191 H], 
[26], [9]. The lambda- lemma says that if finitely many meromorphic functions 
avoid each other in a disc, that is, if their graphs are disjoint, then one can 
always find an additional function which avoids all of them. 

The main conclusion in [T7] is that in any region one can always avoid 
two functions, but in general one cannot avoid three. This is easy to explain 
for the case of avoidance of three rational functions in C. Take the avoided 
functions to be 0, oo and z. If a meromorphic function / avoids them, then it 
must be rational, by Picard's Theorem, but a rational function / that avoids 
and oo in C must be constant, so it cannot avoid z. 

Similar results were obtained in control theory for the case of the unit 
disc or a halfplane. 

We will work in the unit disc from now on. First we give a general 
reformulation of stabilization of three systems in terms of some unusual in- 
terpolation problem. Various special cases of this result are mentioned in the 
control literature, but I could not find a general statement. 

Theorem 1. Let 0i,02 and 03 be three rational functions without common 
poles, and suppose that the set 

E = {z e D : <l>(z) = fc(z) = fa(z)} (4) 

is empty. 

There exists a rational function f which avoids (pi in U if and only if 
there exists a rational function g with the properties: 

(i) divisor of zeros of g coincides with the divisor of zeros of — <p2', 

(ii) divisor of poles of g coincides with the divisor of zeros of </> 3 — <f>\, and 



6 



(in) divisor of ones of g coincides with the divisor of zeros of (pi — (p2- 

Condition that the (pi have no common poles is added only for simplicity of 
formulation: the whole situation is invariant with respect to composition with 
fractional- linear transformations. Condition that E — holds for generic (pi. 

The correspondence between / and g is given by the cross-ratio 

_ (/-0l)(03-0 2 ) 

In the case that there are triple intersections, that is E ^ 0, one has to add 
the condition for each point a e E: 

g{z)(M*) - M*))/(M*) - M*)) = 1 + o(A z ( 5 ) 

where k is the order of the zero of (pi — (p2 at a. 

Thus, simultaneous stabilization of three functions (and the problem of 
avoidance of three functions) is equivalent to finding a function with pre- 
scribed zeros, ones and poles in the unit disc, counting multiplicity, and 
prescribed jets at finitely many points. 

Interestingly, Nevanlinna [22] proposed a similar problem for meromor- 
phic functions in C: to find necessary and sufficient conditions that zeros, 
poles and 1-points of a meromorphic functions must satisfy. Some necessary 
conditions are known [22j [27], see also [251 E3]. Most of these results are for 
meromorphic functions in the plane. 

Consider the following examples. 

1. (Blondel [3]) For which S the following three transfer functions in the 
unit disc are simultaneously stabilizable: 

Pl = z 2 /(z-5), p 2 (z)=z 2 /(z + 5), P3 (z) = ? 

The stabilizer c has to be a rational function without poles in the unit disc 
avoiding l/p%, i = 1,2. This means g(z) = z — c(z)z 2 has to satisfy 

0(0) = 0, g\0) = l, g(z)?±8, \z\ < 1. 

According to a result of Bermant |2J this is possible if and only if 

5 > 5 := 87r 2 /r 4 (l/4), 



7 



and the extremal function is not rational. This inequality gives a necessary 
condition of simultaneous stabilizability. Then an easy approximation argu- 
ment shows that pi,p2 and p 3 are simultaneously stabilizable if and only if 
5 > 87r 2 /r 4 (l/4). 

2. (Patel [21]) For which a > 0, the following three transfer functions in 
the unit disc are simultaneously stabilizable: 

p 1 (z)=z 1 p 2 (z) = z 2 /(z - a), p 3 (z) = ? 

The stabilizer c has to be a rational function without poles in the unit disc 
satisfying 

c(z) 7^ l/z, c(z) 7^ (z — a)/z 2 , \z\ < 1. 
Introducing g = (z — c(z)z 2 )/a we rewrite this in the equivalent form: 

g(z)^l, g(0)=0^z = 0, </(0) = l/a. (6) 

The answer follows from a theorem of Caratheodory (7], [21], [12] 

If a holomorphic function g in the unit disc satisfies (TJJ) then 

\a\ > 1/16. 

There is a real holomorphic function for which equality holds. This extremal 
function is not rational. 

A similar result, but with a smaller constant, was obtained for the first 
time by Hurwitz |15j . 

Now a simple approximation argument shows that the above three sys- 
tems are simultaneously stabilizable if and only if a > 1/16. This answers a 
question stated in [24]. 

Suppose that we wish to stabilize three transfer functions, one of which 
avoids another. The problem is equivalent to finding a rational function 
without zeros and poles in the unit disc, which avoids one rational function 
p. Such c is called a bistable controller. I am not aware of any practical 
application of this "bistability property" by itself, but the desire to control 
three systems with a single controller is reasonable as explained above. The 
problem now is to find necessary and sufficient conditions on a rational func- 
tion p for the existence of c satisfying ()3]) and having no zeros and no poles 
in D. Blondel [3] calls this "one of the major unsolved problems of control 
theory" . 
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A special case of Theorem 1 above, previously established by Blondel, 
says that the problem in equivalent to: 

Finding a rational function w without poles in D, so that 1-points of w 
in D and zeros of w in D are prescribed (with multiplicities). 

We refer to jU [5] and the references in [3] for some necessary conditions 
that zeros and 1-points must satisfy. 

Only one universal restriction (independent of degree) which zeros, poles 
and 1-points of a rational function must satisfy is known. It was found by 
Goldberg pj] and later independently by Blondel. 

To state Goldberg's result, we introduce some notation. Let F be the 
class of all holomorphic functions / in the rings 

p(f) < \z\ < 1, 

with the properties that f(z) £ {0, 1, oo}, and the indices (winding numbers) 
of the curve 

l(f) = {f(z):\z\ = (l + p(f))/2} (7) 

about and 1 are non-zero and distinct. Let F 4 , F 3 , F2, Fx be the subsets 
of Fq which consist of polynomials, rational, holomorphic, and meromorphic 
functions in D, respectively, having finite pairwise distinct numbers of zeros, 
poles and 1-points. We have F 4 C F 2 C F 1 C F and F 2 C F 3 C Fx C F Q . 
The constants p(f) are defined for / G Fj, 1 < j < 4 as 

p(f) = max{|z| : f(z) E {0, 1, 00}}. 

Now we put 

A j = mt{p(f):feF j }, < j < 4. 

Evidently A < Ax < A 3 < A 4 and A < Ax < A 2 < A A . Goldberg's theorem 
says that 

< A = Ax = A 3 < A 2 = A 4 . 

Moreover, extremal functions exist for A and A 2 but do not exist for Ax, A 3 
and At. 

This result shows that if a holomorphic function in the unit disc has finite, 
non-zero, distinct numbers of zeros and 1-points, then these zeros and one 
points cannot lie very close together. 

So we have two absolute constants < A < A 2 which are called Gold- 
berg's constants. The exact value of A Q is known: 

A = exp(-7r 2 /(log(3 + 2y/2)) w 0.003701599, 
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and for A 2 there are estimates 



0.00587465 < A 2 < fi « 0.0252896. 



The constant fi and a function which corresponds to it are conjectured to be 
extremal for A 2 , this function h is described in detail in [6], and we will give 
a short description below. 

If the indices of the curve 7 about and 1 are prescribed to be No, N\, we 
obtain constants Aq(Nq, Ni). One can obtain an exact value of Aq(Nq, N\), 
for any given N x > N > 0, see [5]. 

Being unable to prove that A 2 = fi, the authors of [6] showed that fi is 
the solution of a restricted extremal problem: 

Theorem. A necessary and sufficient condition for the existence of a holo- 
morphic function f in the unit disc, having no poles, a single simple zero at 
a and a single multiple 1 -point at —a is that \a\ > fi. If a = fi this function 
is unique and transcendental. If \a\ > fi there exists a polynomial f with the 
stated properties. 

Thus in the simplest case of one simple zero, one multiple 1-point and 
no poles, we have a necessary and sufficient condition for the existence of a 
rational function with prescribed zeros, 1-points and poles in the unit disc. 

This can be restated as a necessary and sufficient condition for a stabi- 
lization problem as follows: 

The transfer function 



can be stabilized by a bistable controller if and only if \a\ > fi. 

In more complicated cases, there is no hope for such simple conditions. 
For example, Blondel [3] states the following problenjl: 

For which 5 > there exists a rational function which in the unit disc has 
no poles, a single simple zero at 0, and exactly two simple 1-points ±i5? 

It is known that there exists So > 0, with the property that such function 
exists for 5 > 5 and does not exist for 5 < 5 . 

2 He even offered a prize of 1 kg of fine Belgian chocolate for this problem. Nevertheless 
it is still wide open. 




z — a 
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Evidently 5q > A 2 and it is not difficult to show that this inequality is 
strict. The current world record [5] for the estimate from above seems to be 
5q < 0.1148. The best known lower estimate is 0.01450779. It can be obtained 
from the estimate in jTJ] of the minimal length of a closed hyperbolic geodesic 
in a twice punctured disc [6]. 

In conclusion, we sketch the definition of the function which is conjectured 
to be extremal for A 2 . The fundamental group T of C\{0, 1} is a free group 
generated by simple loops A and B around and 1. Let V be the subgroup 
generated by A and B 2 . It is also a free group on two generators. Let 
g : X — > C\{0,1} be the covering map corresponding to this subgroup V, so 
that T' is the fundamental group of X. One can show that X is a Riemann 
surface which is conformally equivalent to the twice punctured disc, and 
we can identify it with D\{— fi, fi} for some /x G D. Then g becomes a 
holomorphic function in D which has one simple zero, say at — fi and one 
double 1-point at ji. We conjecture that A 2 = \i. One can express our function 
g in terms of solutions of a Lame equation and modular functions. 
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